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We present a formalism to determine the imaginary part of a general chiral model 
in the derivative expansion. Our formalism is based on the worldline path integral 



for the covariant current that can be given in an explicit chiral and gauge covariant 
form. The effective action is then obtained by integrating the covariant current, 
taking account of the anomaly. 
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■ I. INTRODUCTION 

o 

>' 

When discussing the influence of fermions on the dynamics of some theory, e.g. of the 
Standard Particle Model (SM) in a cosmology setting, it is mandatory to integrate out the 
fermions. This procedure in onedoop order already becomes quite involved if arbitrary chiral 



couplings o 
of Salcedo 



: sp_ace-time dependent outer bosonic fields are considered. The extensive work 
1 



resulted in an effective action to leading order in covariant derivatives of 
such fields. Basis is a refined calculation in momentum space, handling of anomalies using 
the Wess-Zumino-Witten model (WZW), and last but not least, a very practical shorthand 
notation. 
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Such effective actions are quite important in evaluating some models. For example, in 
ref. B Smit discussed a form of 'cold' electroweak baryogenesis at the end of electroweak 



scale inflation 



4j which could very well work if the rephasing invariant 
J = s 2 1 s 2 s 3 c 1 C2C 3 sin(<J) = (3.0 ± 0.3) x 10~ 5 



(1) 



of the Jarlskog determinant is not accompanied by further suppressions through mass ratios. 
It was proposed 3_| that derivative terms in the effective action that are analytic in the time- 
dependent masses considered non-perturbatively could be very important in non-equilibrium. 
Such effects were also observed in ref. In the work 3] the fourth order derivative result 
of ref. JJ turned out not to contain CP violation, but the claim was made that higher orders 
of the imaginary part of the effective action will do. 

Worldline methods in first quantized quantum field theory are ideally adapted for cal- 
culating effective actions: One considers the propagation of a particle in some space-time 

n aa, 



dependent background [6] , but in x-space path integral formu 



also related to the infinite tension limit of String theory 
discussion of various effective actions in one-loop [9|. 
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ation [7|. This method 



was used heavily for the 
14| and two-loop 15, 1J, Q, Q 
order. For example, the high order in the inverse mass calculation of ref. [15J could hardly 
be done with other methods. 

The present paper provides a formalism to determine higher order contributions to the 
imaginary part of the effective action using the worldline formalism. We are concerned with 
the effective action of a multiplet of N Dirac fermions coupled to an arbitrary matrix- valued 
set of fields, including a scalar $, a pseudoscalar II, a vector A, a pseudovector B, and an 
antisymmetric tensor K^ v . One peculiar feature of the imaginary part of the effective action 
is that it cannot be written in a manifest chiral covariant way, due to the presence of the 
chiral anomaly. One possibility to arrive at a closed expression for the effective action is to 
abandon manifest chiral covariance as it was done in ref. The resulting expression is 
rather complicated and not well suited for higher order calculations. Alternatively, it was 
proposed in ref. [lj to determine the covariant current for which a manifestly chiral covariant 
expression exists and to take account of the anomaly when integrating the current to yield 
the effective action. Following this idea, we present a worldline path integral formulation of 
the covariant current. 

Before we do so, we review the worldline formalism by discussing the derivation of the 
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real part of the effective action. A single Dirac fermion in the presence of both a scalar and 



pseudoscalar fields in the context of the worldline formalism was first treated in ref. 



20|, 



the inclusion of a pseudovector in ref. 2l(] . In our discussion of the real part of the effective 



action we will follow the elegant subsequent work of refs. [19|, l22j ]. 

Section [III contains the derivation of the real part of the effective action, the derivation 
of the covariant current and the matching procedure to obtain the imaginary part of the 
effective action. In section HTT| we briefly reproduce the results in lowest order from ref. 
As a novel result we present the imaginary part of the effective action in two dimensions in 
next to leading order in section [IVl 

II. EFFECTIVE ACTION 

We are concerned with the effective action 

iW[&,n,A,B,K\ = logDeti[i#-$ + i 7 5 IT + ^ + 7 5 $ + ij^K^], (2) 

and its continuation to Euclidean space. The 7 matrices remain unaffected by the contin- 
uation, but it is useful to introduce the following notation, (je)] = ijj, (7^)4 = 7o ; an d 
(7^)5 = 75- After Wick-rotation, t —>■ —it, one obtains with this new notation 

-> i0 E , 4 -> Ue, $ -> i$ E , -f-fK^ -> -{7e)M„K e ^. (3) 

From now on, the E subscript will be suppressed. 
The effective action of Eq. (j5J) now reads 

-W[^,U,A,B,K]=logBet[0], (4) 

with the operator O in momentum space defined by 

O = y - i$(ar) - 75 n(x) - £(x) - 7 5^(x) + J^K^. (5) 

As in ref. {2], Q], the real and imaginary parts of the effective action are analyzed separately 

-W + -iW = log(|Det[C]|) +i arg(Det[C]). (6) 



A perturbative expansion in weak fields [22| shows that graphs with an even number of 75 
vertices are real, and graphs with an odd number of 75 vertices are imaginary. This will 
prove useful when the behavior of the effective action under complex conjugation is explored 
later on. 
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A. Real Part of the Effective Action 

Our intention is to obtain a worldline representation for the effective action with manifest 
chiral and gauge invariance. This is unproblematic for the real part, but it causes certain 
difficulties for the imaginary part due to the chiral anomaly. In order to familiarize the reader 
with the worldline method we review the derivation for the real part in four dimensions as 



it was presented in refs. 
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1. Construction of a Positive Operator for the Real Part of the Effective Action 

In order to use the worldline formalism, one has to rewrite the effective action in terms 
of a positive operator, thus obtaining 

W + = -ilogDetffltO]. (7) 

The problem with this operator is that it contains terms linear in the 7 matrices, what makes 
the transition to a path integral of Grassman fields problematic. One way to avoid this 
problem is by doubling the fermion system and exchanging the operator O for a Hermitian 
operator £ yielding 

^+ = _Il ogD et[OtO] = -ilogDet[S 2 ], £ = f ° ° I . (8) 

2 4 \C?t J 

Since £ is Hermitian, one can use the Schwinger integral representation of the logarithm 
without any restrictions. One obtains 



W 



1 f°° rIT 

lJ o ^Tr exp(-TS 2 ). (9) 



At this point, it is natural to introduce six 8x8 Hermitian Ta matrices. These matrices 
satisfy {Ta, Tb} = 2o~ab, with A, B = 1..6 and are defined as 

r„= I" r 5 = I 75 1 , r 6 = I " -|. (10) 




v 7m J yy 5 

For later use we also introduce the equivalent of 75, 





5 



and T 7 anticommutes with all other T matrices. 

Expressing E in terms of these new matrices yields Q, 

e = r> M - A^) - r 6 $ - r 5 n - % r M r 5 r 6J B M - % r fl r u r fi K t 



hp • 



(12) 



The aim is to turn Eq. ffl~2l) into an expression which is manifestly chiral covariant. This 
can be achieved by changing to a basis in which i T^Tq is diagonal 22] using the following 
transformation 




M 



fh o o o\ 

h 
h 

\o h o 0/ 



(13) 



In this basis, £ takes the form 



>; = M _1 SM 



. T5 (_ ijff t + i TA4Tj/ ^-t 



-iH + \ llllv Kl 
ln(Pn ~ A*) 



which is manifestly chiral covariant. Here A L = A + B, A R = A — B, H 



(14) 
$ - iU, 



K s = K - i K and K 



-^t^ LUp(7 K fXJ have been defined. 



The square of £ constitutes a positive operator which is suitable for the worldline formal- 
ism. However, even though this expression contains only even combinations of 7 matrices, 
the coherent state formalism cannot yet be used to transform this expression into a fermionic 
path integral. In the coherent state formalism, the 75 matrices have to be rewritten as a 
product of the other 7 matrices, what would result again in odd combinations. One possi- 
ble solution of this problem is to enlarge the Clifford space, replacing the 7 matrices by V 
matrices 

x o 1 

1 

The matrix T 5 is then independent from the other T matrices and the coherent state for- 
malism with six (instead of four) operators can be used. The doubling of the Clifford space 
inside the trace has to be compensated by a factor \, such that Eq. ([9]) reads 

1 



1A 



1A 



A E [1 ... 5]. 



(15) 



W 



^-Tr exp(-TS 2 ) 



(16) 
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and the operator £ 2 is given by 

t 2 = (p - A) 2 + h 2 + \k^ v + ^r„r„(.7> + {H, JC^} + % [K m K pv \) 

+i Ff.TsCDf/H + {p u - A u , /C^}) - -r (U/0(T r 5 P M /C p(T - -T^pvtC^tCpv, (17) 

with enlarged background fields defined by 

I A a °\ / iH\ I iAT* , 

r^...^ = r^.-.r^] denotes the anti-symmetrized product of k Y matrices, and the field- 
strength and the covariant derivative have been defined as 

J> = dpAv - d u A^ - i [A,,, A v ], V^x = d fl x~i [\, x}- (19) 

The S 2 operator is seen to be manifestly gauge and chiral covariant. It also contains T 
matrices to even powers only, and is well suited for the worldline path integral representation. 



2. Worldline Path Integral 



With the use of the coherent state formalism 22|, |23j, one can perform the transition 
from T matrices to a path integral over Grassman fields if), with the correspondence T^r^ — > 
Zip Aip b an d TaT b T c T£) — > AipAipBipcipD, as long as A, B, C, and D are all different. The 
final form for the real part of the effective action is 

W + = - I™ *^-N I Vx I V^tiVe-^ drC{T \ (20) 
8 Jo T J J AP 

Here M denotes a normalization constant coming from a momentum integration and AP 
stands for antiperiodic boundary conditions, which must be fulfilled by the Grassman vari- 
ables ip{T) = —ip(0). The Lagrangian is given by 

x 2 1 1 
£(t) = — + -ipAi>A ~ i XpAp + U 2 + -JC^JC^ + 2i 4>^ 5 (V^H + i x u IC^ u ) 

+ i [T^ + {H, /C^} + i [/C w , K pv \) 

- ^^Ipvilp (2^5^/C^ + IpafcfivKpa) • (21) 



The periodic boundary conditions for the field x{r) suggest to separate the zero modes of 

(It 2 ■ 



the free field operator The fields x(t) are split into a constant part and a r dependent 
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part according to x(r) = x + y(r), with d T x = and J Q dry{r) = 0, and the measure in 
the integral is changed into Dx = Dyd D x§. The Green function is defined on a subspace 
orthogonal to the zero modes. The ipA fields contain no zero modes, so that the propagators 
for the y{r) and iPa(t) fields read 

(y(ri)y(r 2 )) = ^— ^ -\t x -t 2 \, 
(^ACrOVtefe)) = ^Assign (n - r 2 ) . (22) 

This formalism can then be used to determine the real part of the effective action as discussed 
in ref. M. 



B. Imaginary Part of the Effective Action 

As in the case of the real part of the effective action, one requires a positive operator in 
order to use the Schwinger trick. Even though this is still possible for the imaginary part, 



gauge and chira 
example, in ref. 



invariance cannot be manifestly conserved due to the chiral anomaly. For 



19 



22J a parameter a is introduced, which breaks the chiral invariance, but 
leads to a positive operator. However the resulting expression is not appropriate for higher 
order calculations since the breaking of manifest chiral invariance leads to a large number 
of contributions in the perturbative expansion of the path integral. 

The aim of the present work is to present a worldline representation of the effective 
current for which a manifestly chiral covariant expression exists. This current can then be 
integrated to obtain the effective action [l|, Q, 25]. This integration rather proceeds by 



matching: First, a general effective action is proposed, which has the expected chiral and 
covariant properties. The functional variation of this action is then matched to the covariant 
current that is obtained using the worldline formalism. This method has the advantage that 
it is both gauge and chiral invariant at each stage of the calculation. The anomaly only leads 
to additional complications in the matching procedure of the lowest order contributions as 
will be discussed in detail in the next section. 

Starting point of our analysis is the functional derivative of the imaginary part of the 
effective action in Eq. (jHJ) 

SW~ = ^5 (log Det O - log Det C f ) = ^Tr (tiO^ - SO^j . (23) 
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This expression can be rewritten in terms of a positive operator which can be used to employ 
the worldline representation in combination with the heat kernel formula. Incidentally, it 
can also be expressed in a manifestly chiral covariant form, what simplifies higher order 
calculations tremendously as compared to the formalism presented in ref. [3]. 



1. Construction of a Positive Operator for the Imaginary Part of the Effective Action 



The expression in Eq. (1231) can be transformed using the operator E defined in Eq. ([8]) 

1 ( 50\ ( 1/cA 

SW~ = -TV 1 , (24) 

2 \-5<y J \l/0 J 
which, with the introduction of a new matrix x, can be rewritten as 

5W- = ^Tr X SZZ-\ (25) 

with 

' o o\ /mo, 

ct o / yo -h 

To produce the positive definite operator E 2 in Eq. (1251) . we multiply and divide by E, using 
the cyclic property of the trace and the fact that E anticommutes with x, to obtain 

SW~ = ^Tr (x^SE + Ex5E) E" 2 

= ^Tr X [5E,E]E' 2 . (27) 

Since the last factor is a positive operator, it can be reexpressed as an integral, similar to 
the expression of the real part of the effective action in Eq. (fT7|) . namely 



i r°° 

SW~ = -Tr / cffx[5E,E]e- Ts2 . (28) 
4 Jo 

As in the case for the real part, the chiral covariance can be made manifest by changing to 
an appropriate basis. With the help of the matrix M in Eq. (TlBl) . one obtains again 

£ = 7 M (JV ~ A) ~ 75^ - ^7/^75^- (29) 
The additional factors x [<^, E] read 

M- 1 X M = x= f 75 ° | =X7s, (30) 
-7s. 
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and for the case ST, = —J^SA^ 

+i 757/. [8A V , Kfu,] ~ 2 757^a ct {SA^, lC X a} ■ (31) 

To use the coherent state formalism, it is again necessary to enlarge the Clifford algebra 
and to replace the 7 matrices by T matrices. However, taking into account the factor 75 
in Eq. (1301) the imaginary part of the effective action contains only odd combinations of 7 
matrices. Thus, the replacement 




1a -> T A = lA ® I I , A e [1 ... 5] (32) 
has to be compensated by a factor 



1 



i 



- 2 r rr 6 = 14 ® I 1 q I ■ ( 33 ) 
The overall factor T 7 changes the boundary condition of the fermionic sector from antiperi 



odic to periodic as explained in ref. 22J. This means that the fermionic sector contains zero 
modes, which have to be separated in the same way as was done for the bosonic sector. 

Including the factor in Eq. ( 1331) to compensate for the doubling of the Clifford space, one 
obtains 

5W~ = -Tt dTT 7 T 6X w(T)e- TJ]2 } (34) 
where S 2 is given in Eq. (jTTI) . and the insertion due to the commutator yields 

w(T) = —TsT^iSA^Pv-A^-iTsV^A^-T^SA^H} 

+i [5 A V1 JC^} - -T^a {$A^, IC X a} ■ (35) 

To transform this expression into a worldline path integral, a similar procedure as for the 
real part of the effective action can be followed. Products of T matrices can be replaced by 
Grassman fields, however in this case the Jacobian of the transformation contains additional 
contributions from the zero modes 

veve = de^dOxdhdhdhve've' 

= -#i#2#3#4#5#6^'- (36) 

'J 
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The factor J only includes the Jacobian for the zero modes, while the Jacobian for the 
orthogonal modes is absorbed in the normalization of the correlation functions of the ip' A . 
J can be calculated from the definition of the Grassman fields ip in the coherent state 



J = deb l^f) =H)( B + 2 V 2 . (37) 



formalism 22j and yields in D dimension 

9 

dip 

The final result can be expressed as 

5W~ = ^tr J dTAf Jvxf T>ipX™(T)Ve- f ° dT£( - T \ (38) 

The Lagrangian is of the same form as in the real part, Eq. (I2TT) . 

x 2 1 1 
C(t) = — + -vpAipA - i XfiAft + H 2 - -K^Kftu + 2i ip^s (V^H + i ± V K.^) 

and the trivial integration over ip§ can been carried out, so that the insertion yields 

w(T) = -Ai 4> 5 4>^ u 5A^x u - 2i faV^SA^ - 2^ {5A^, H} 

+2i [5A V , JC^} - 2i {SA„, /C Aff } . (40) 

The normalization Af coming from the momentum integration, satisfies 

AT J Vxe~£ dT i = (4vrT)~ D/2 J d D x. (41) 

The Green function for the bosonic field x is the same as for the real part of the effective 
action, Eq. ( 1221) . while the Green function of the Grassman fields ipA differs due to the 
presence of the zero modes. The fermionic fields are split according to iPa( t ) — ^ + 1 I j 'a( t )^ 
with d T ip A = and J Q T drip' A (r) = and the measure turns into T>ip = dipid^d^d^d^'D^j' . 
The Green function for the ip A fields, defined on a space orthogonal to the zero modes, reads 

(^(n)^(r 2 )} = S AB QsignCn - r 2 ) - ^j^j . (42) 

These results can be easily generalized to different dimensions. In two dimension, one obtains 
an additional overall factor — i from the Jacobian of the zero modes and the fermionic 
measure reads Vip = dvpid^dvpsVip' . 
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2. The Effective Density 

The effective density is obtained by varying with respect to the TC field, so that <5X = 
— 75 In comparison to the worldline representation of the covariant current only the 
insertion changes into 



7s7m {SH, p lt -A lt }+ [SH, H] + % - lillv [SH, K^] . (43) 



The corresponding insertion w (T) in the path integral reads then 

w(T) = -2% ^x^SH + 2^5 [SH, H] + 2i ^if u [SH, . (44) 

Since 5 A carries an index, the effective current is of one order lower than the effective density 
and usually results in less terms to calculate. The advantage of the effective density lies in 
the matching process, since the factors in the effective density consist of the same type as 
found in the effective action. They both combine the same type of object, T>H and T , to the 
same kind of order, while the effective current combines the terms to a lower order. Besides, 
there is no distinction between a consistent effective density and a covariant effective density, 
as there is for the effective current, as will be explained in the next section. 



3. Distinction between the Consistent and the Covariant Current 

With Eq. fl38l) an expression for the covariant current which is chiral and gauge covariant 
was derived. This current cannot be the variation of the effective action, since the effective 
action contains the chiral anomaly, and in fact the covariant current is not a variation of 
any action. The reason for this is that performing the variation does not commute with the 
regularization procedure we used, namely the Schwinger trick. On the other hand, knowing 
the chiral anomaly, one can reproduce the so-called consistent current that denotes the true 
variation of the effective action. 

To explain the relation between the two currents, we define a general variation 

s - = h Y >^y (45 > 

so that a gauge variation 5^ is given by 
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Two subsequent variations have then the commutator [5y,<y = S[y,€\ an d m order to find 
the transformation properties of the consistent current, one can apply this commutator to 
the effective action 

[5 Y , S^W-iA,] = 5 [Y ^W-[A,]. (47) 



Using the anomalous Ward identity [26|] 

S^W-lAp] = J dx^(x)G[A^](x), (48) 

with G[^4 M ] denoting the consistent anomaly, one can evaluate both sides of Eq. ( f471) to 
obtain 

J dx[Y^i}(x)jJ-^W-[A,} = 5 Y J dxax)G[A,]{x) 

-6 e J dxY,{x)j^W-[A,]. (49) 
Defining the consistent current as the variation of the effective action 

= im w ' [A - ] ' (50) 

one finds 

J dxY^x)5^{f{x)) = J dxY„[(f{x)),t](x) + J dxi{x)5 Y G[A„]{x). (51) 
Since Y was a general variation this leads to 

WW) = Wix)),*] + f d y^y)jJ^) G lA)(y)- (52) 

This shows that only if the anomaly vanishes, the current transforms covariantly. This 
relation can be used to determine the connection between the consistent current, i.e. the 
true variation of the action, and the covariant current. The latter is obtained by adding an 



object P^f-A^], called the Bardeen-Zumino polynomial 271 ]. to the consistent current so that 
the sum transforms covariantly 

m = (f) + (p»)- (53) 

This implies the following gauge transformation property for the BZ polynomial 

8^[A,\{x) = [PI^U] (*) - J dyay)jJj^G[A,}(y). (54) 



13 



It is not obvious that such an object exists, but using 



1 



48tp 



e^ti X {A v Txa + F\oA v + % A v A\Aa, 



and the consistent anomaly 26] 



G[A, 



1 



e^ Xa tT X da I AvdxA* - ^A v A\A a 



(55) 



(56) 



it can be shown that the definition of P M in Eq. ( 1551) provides a unique polynomial in A 
that satisfies Eq. (1541) . The corresponding functions in two dimensions are given by 

P" = ^-e^tr A, G[AJ = -^e^tr X A- (57) 

47T 47T 

As stated above, the path integral in Eq. ( |38l) constitutes a worldline representation of 
the covariant current. To obtain the imaginary part of the effective action from the covariant 
current one can use the following ansatz 



gWZW 



+ w: 



(58) 



Here, T g wzw is an extended gauged Wess-Zumino-Witten action [24], |28|, |29( , which is chosen 
to reproduce the correct chiral anomaly, and W~ denotes a chiral invariant part. The 
variation of the functional T g wzw, consists of a part that saturates the anomaly, namely the 
BZ polynomial, and a covariant remainder which has to be added to the variation of W~ to 
yield the covariant current. 



4- The Wess-Zumino-Witten action 

When the effective action is separated into two parts, it is required by the non-covariant 
part that it reproduces the anomaly. It is well known that the WZW action has this property. 
The ungauged WZW action in four dimension is e.g. of the form 



T(U) = — -J / d 5 xe abcde tr 



Q 



5 



(59) 



48vr 2 

where Q is a five-dimensional space with boundary dQ equal to the R 4 flat Euclidean 
space. The matrix U is a unitary matrix, and is usually related to the case where the 
mass can be expressed as a constant times that unitary matrix. We are interested in the 
more general case when the mass matrix is not of this form which is called extended WZW 
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action. In addition, the presence of the background gauge fields makes a gauging of the 
action mandatory. The gauged extended WZW action can be generally expressed as the 
integral in five dimensions 24] . Unlike the action itself, the resulting current turns out to 
be a total derivative in five dimensions, such that it can be represented by an integral over 
the physical four- dimensional space 

5T gWZ w = 9^5/ d i xe^ Xa tT X [SA fl (-n- 1 V v H'H- 1 V x 'Hn- 1 V (T 'H 

+v u nn- 1 v x nn- 1 v a nn- 1 - % {n- l v v n - v v nn- l ,T Xa ] 

+ l -H{H- 1 V u H,F X(7 }H- 1 - l -U- x {V v HH-\T x *}n 

-2 {A, Fx*} - 2i AuAxAa)] , (60) 

or in two dimensions 

STgWzw = / d 2 x e^tr X [SA^ i-i H^VJH + % VJUH^ - 2A V )\ . (61) 

Notice that in both cases the last term in the current denotes the BZ polynomial. The 
remaining chiral covariant terms have to be subtracted from the covariant current before it 
is matched to the effective action according to the ansatz made in Eq. (!58lh 



III. LOWEST ORDER EFFECTIVE ACTION 



A. Effective covariant current 

In order to reproduce the results from ref . [l| , we neglect in this section the antisymmetric 
field K^y. The fields A and Ti are matrices of some internal group, and we only assume that 
7i,(xo) is nowhere singular. With this in mind, we restate our result Eq. ( 1381) from the last 
section in D dimensions 

•D/2 poo /• /• T 

5W~ = — — tr J dTAf J Vx J Thp xwfflPe~ So drC{r \ (62) 

with 

x 2 1 

C(t) = — + -iPa^a - i XpAp + H 2 + 2i ^^V^H + % i^^J 7 ^, 

w(T) = -Ufa^JApX v -2i$ B Vp8A ll -2%{5A llt H}. (63) 

Next, the derivative expansion of the heat kernel is used. In the derivative expansion terms 
are classified by the number of covariant indices that they carry, so that T>^7i is of first order, 
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while is of second order. The worldline formalism is well suited for this expansion, and 
there are two major advantages compared to the more traditional methods used e.g. in 
ref. First, the tedious manipulations using the 7 algebra are avoided. Secondly, the 
momentum integration is omitted and replaced by the rather trivial integration in r space. 



The coordinate is split as x(t) = xo+y(r), and we work in the Fock-Schwinger gauge 30j, 
in which A(x) ■ y = 0. In this gauge, expressions remain gauge covariant and the field A can 
be expressed in terms of the field strength tensor T pv by 

A ii {x)= \ daaJ r pp (x + ay)y p . (64) 
J 

All background fields can then be expanded around the point xq in terms of covariant 
derivatives 

X{x + y(r)) = exp (y(r) ■ V xo ) X(x ), (65) 



where T> XQ refers to the covariant derivative in Eq. ([19]) with respect to x . With the 
expansion of the field strength tensor in terms of covariant derivatives and Eq. fl64l) . one can 
rewrite the field A as 

Api?) = ^y P ^( x o) + ^VaVp'DaFpiiM + -^^y a ypy p V a V p F p ^x ) + ... . (66) 

Since we will not carry out the integration with respect to xq we use the following notation 
in D dimensions 



D/2 



(i) trx / dDx ° x - (67) 



It is important to remember that x an d Ti. anticommute; hence, when the cyclic property of 
the trace is used, a minus sign is generated, for example 

= -(e"^ 3 ^}. (68) 

After expanding the mass field Ti(x) 2 = H 2 (xo) + y p V p 7{ 1 {xQ) + . . ., the field H(xo) is 
treated non-perturbatively. Since all the fields can be matrices of some internal space the 
resulting expressions normally cannot be expressed in closed form. For this case we use the 
labeled operator notation laid down in ref. 0, [y]]. The notation works as follows: In an 
expression f(Ai, B 2 , . . .)XY . . ., the labels of the operators A, B, . . . denote the position of 
that operator with respect to the remaining operators XY .... For instance, for the function 
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f{A,B) = a(A)/3(B), the expression f(A 1 ,B 2 )XY represents a(A)XP(B)Y. In the case 
at hand, the operator appearing in the functions is always m := TI(xq), such that general 
functions / can be easily interpreted in the basis where m is diagonal. Using this notation, 
Eq. ( 1681) can be recast as 

(e^HF, u H 3 F Xa ) = (f^m^lT^T^) = - (e^ Xa m 3 ml^ Xa ) 

= _ ( e ^m 2 ml^ Xa ) = - (e^H 3 F, u HF Xa ) . (69) 

This notation can also be used to simplify the matrix valued derivative. Using the definition 

fT7f \f \ /K) - /(m 2 ) , . 

(V/)(mi,m 2 ) := , 70 

mi — m 2 

it is possible to prove that 

P M /(m) = (V/)(mi,m 2 )XV,«. (71) 
For example, in the polynomial case f(m) = m 3 one obtains 



VJ(m) = V ll (H i )=V^HH 2 + HV ll HH + H 2 V tl H 

ml — m\ 

m i — m 2 

(Vf)(m 1 ,m 2 )V^H. (72) 



(ml + m x m 2 + m?) V^H = ^ — ^ V^H 



As mentioned earlier, non-polynomial expressions are hereby interpreted in a basis where m 
is diagonal, so that for m = diag(di, . . . , d n ) 

f(mi) - f( m 2) x = f(dj) - f(dj) x 
m i — m 2 d-i — dj 13 

More general, this suggests the following definition for the case with several variables: 

v7 ,/ x f (m u ...,m k+u ...,m n ) - f (m u ...,m k ,...,m n ) 

Vfc/(mi, . . .,m n ) = , (74) 

m k - m k+1 

where m k indicates that the corresponding argument is left out. 

If all arguments of the functions are of the same type one can further simplify the notation 
and use subscripts to refer to the argument of the function, e.g. f(mi,m 2 ) =: / 12 and we 
employ this notation in the following. Additionally, negative arguments will be denoted by 
underlining the corresponding index, /(— 7711, m 2 ) —'■ fi2- More applications of the labeled 



operator notation can be found in refs. [l 
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The path ordering in Eq. (1621) is defined by 

N pT pT pt\ ^tjv-i pT pT _~} ~ 

V\\ dn = N\ dr x dr 2 --- dr N = N\ dr x ■ ■ ■ / dr N TT Ofc - r i+1 ). 

i=l JO JO Jo Jo Jo Jo 

(75) 

Separating the Lagrangian Eq. ( 163]) into £(r) = £o( r ) + 7^ 2 (^o) + £i( r )j with 

x 2 1 
£o{t) = — + -ipAipA, 

Ci(t) = -i ipAp(x) + 2i ip^V^Hix) + i ijj^T^x) + y IJ T> l Ji' l {x Q ) + (76) 

the terms of the expansion of H 2 (x), except the leading term H 2 (xo), are attributed to £i(t), 
and treated perturbatively. Notice that Cq commutes with the rest of the Lagrangian, so 
that the expansion of the path ordered exponential in Eq. (162|) takes the form 

+ dn jH dr 2 e^ T -^ HH ^ (-£i(n)) e -(— 2 )« 2 (-o) (-^(t-2)) e-^w 2 ^) + . . ^ . (77) 

When performing the ip integrals, the zero modes have to be saturated and at least a factor 
t/)° . . . ?pE>?p5 is required from the Grassman fields in order to contribute. The first term in 
Eq. (ITT)) lacks the appropriate r/> factor except in two dimensions, where the first term of 
the insertion Eq. (163]) already has the appropriate factor. However it contains a factor x^ 
which must be contracted with a similar factor to form a Green function, hence it does not 
contribute and can be left out. The rest of Eq. ( 1771) can be simplified using the labeled 
operator notation. Using the expression rr? n to denote TC 2 (xo) in nth position, one obtains 

Ve -f T drC(T) = e -J T drC (r) f ^ ^-^(^j^^ 

+ £ dr i £ dne-^M^-^M^-^Ctin)^^) + ...V (78) 

The evaluation of the worldline path integral can be summarized as follows: First, all fields 
in Eq. (1751) and the insertion are expanded around Xq. Next, the functional integration over 
the y fields is carried out, generating bosonic Green functions. Then, the ip integrations are 
performed saturating the zero modes and generating fermionic Green functions. Finally, the 
T and r integrations are performed. 

Before presenting the actual calculation, we comment on the behavior of the effective 
action under complex conjugation. As noted earlier, in any contribution to the imaginary 
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part of the action the field ip^ appears an odd number of times. If one attributes a factor 
i to the operators T and 5A one observes that the remaining expressions in the current 
in Eq. (1631) are real. Accordingly, all expressions in W~ are real as long as a factor i is 
attributed to the operator T . In addition, notice that the effective action has to be an even 
function in the masses due to chiral invariance. 

In order to showcase the method, we present the lowest order calculation in two dimen- 
sions. The lowest order contribution coming from the first term in the insertion is given 

b y T 

-Ai^^ u y v (T) [ dne- Tm *-^-™l)y a ( T )V a H 2 5A,(T). (79) 
Jo 

Performing the y and i/j integrals one obtains 



m 1+ m 2 ) — dne-^M^-^gBp^VpHSA* 
Jo 1 Jo 

4" Jl 2 {m x + m 2 ) V^H 5A V ) . (80) 



i 

~ 2 

The second term of the insertion does not contribute at lowest order since it is already of 
second order in derivatives but lacks the appropriate fermionic factor to saturate the zero 
modes. The third term of the insertion leads only to one contribution of the form 

-2%{5A^H} f dne- Tm *-<™%-™l) {-2i^ b V u H). (81) 
Jo 



yielding 



2 \ Jo T J 

= ~\ (f v JUrnx ~ m 2 )V,H5A v ) . (82) 

The factor (m x — m 2 ) results from the anticommutator in Eq. (1811) . and the sign change 
in the cyclic property of the trace as explained in Eq. ( 168|) . The integrals J are given in 
Appendix [A] The total current is hence given by 



5W~ = -i ( e^A\ 2 V„U 5A V ) , (83) 

A i = 1 m!m 2 log (ml/ml) 

m\ — vri2 (mi — ?7i 2 )(mf — m^) ' 



where the function A\ 2 has been defined. This agrees with the results obtained in ref. 
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B. Effective Density 

The effective density can be obtained similar as the covariant current, utilizing the inser- 
tion in Eq. (liljl . Neglecting the antisymmetric tensor JC, the insertion is 

w(T) = -2i 4>^x^H + 2^ 5 [5H, U\ . (85) 

The contributions to the effective density are 



~\ (^123( m l + m 3) + ^123( m l + m 2) - ^l 7 23( m 2 + m 3 )) V^HVJtj SH 

^ (y-\B\ 2 F^ + B 2 123 V^HV u 7?j 6H^ . (86) 



where the functions _Bi 2 and -B123 are given by 



with 



$> = ; ^ log(i), (87) 

rrii + m 2 (mi + m 2 )(mf — m 2 ) \ m 2/ 

B 2 123 = B« 3 + B[ 23 log(m?) + Bg, log(m 2 2 ) + Bg 2 log(mj) , (88) 



^123 = 7 w w 7; (89) 

(mi - m 2 ){m 2 - m 3 J(mi + m 3 ) 

b l = (mf + mim 2 m 3 ) 

12 |L (mi — m 2 )(mi + m 3 )(m^ — m^mj + mf) ' 



in accordance with ref. |l|]. 



C. Effective Action 



We proceed and briefly present the derivation of the imaginary part of the effective action 
following ref. pj. Using the ansatz in Eq. (!58j) . the most general functional for W~ consistent 
with chiral and gauge invariance in two dimensions reads 

W~ = (frNnVpKDv'H) . (91) 

An additional term proportional to T could be added but it can be removed by partial 
integration. Notice that N12 is a real function according to the comments made in the last 
section. 
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The function Ni 2 has some nontrivial restrictions. First of all, the function Ni 2 is even 
in m such that 

N(-mi, -m 2 ) := Nu = N 12 . (92) 
Because of the cyclic property of the trace one obtains 

N 12 = N32 = N g i = N21, (93) 

and due to the Hermiticity of W~ 

N 12 = -N 32 = -N l2 = -N 2X . (94) 

Varying VF c - [.4, Tt] with respect to A, one obtains 

5W~ = -1 (? v (-2 (mi + m 2 )N l2 ) V^H5A U ) . (95) 

Comparing this to Eq. (1531 and adding the covariant contribution in Eq. flBTj) coming from 
Fgwzw one has 

1 mim 2 \og(m\/m 2 l ) 1 1 



vn\ — m 2 {nil — m 2 )(m\ — m^) 2rrii 2m 2 
which finally leads to 



2 (mi + m 2 )N l2 , (96) 



N^ = \^^ 2 ( X ^^-\{\ + - 2 )). (97) 

lm\ — 777-2 \ m l ~ m 2 ^ \mf 777^ / / 

At higher order, the matching of the effective potential to the current potentially becomes 
more intricate. On the other hand, the anomaly only contributes to the leading order, 
such that the knowledge of the covariant current (that in higher order coincides with the 
consistent current) suffices to determine the effective action. 



D. Four Dimensions 

For completeness, we also present the results for the effective action and the effective 
current in four dimensions. The matching procedure proceeds the same way as in ref. 
and we do not repeat it here. 

The effective current in four dimensions consists out of three terms and reads 



(98) 
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while the effective density can be written as 

1 % 

- Bj 334B v IM nv v nv x nv <r 'Hj sn^ . (99) 

The functions A\ 23 , A 3 123 , Af 2U , Bf 23 , Bf 2U , Bf 2U , Bf 234 , and B 7 12345 are given in Appendix[Bl 

IV. NEXT TO LEADING ORDER EFFECTIVE ACTION IN TWO 

DIMENSIONS 

In this section we present as a novel result the imaginary part of the effective action in 
next to leading order and two dimensions. Even though the results are rather lengthy, the 
evaluation of the worldline path integral involves only very basic integrals such that it can 
be easily implemented using a computer algebra system. 

In two dimensions and in next to leading order, the imaginary part of the effective action 
takes the form 

+Ri23V a V a HV fJ/ HV u H + ^R\2zT ^ V V oiHV oiH 
+M 1234 V f ,HV a HV u HV a H)} . (100) 

At next to leading order the action is chiral invariant and the effective action can hence be 
immediately obtained by matching with the covariant current that in this order coincides 
with the consistent current. These functions must have the following properties 

Pl2 = -Pl2=P21, Ql2 =Ql2 = -Q21, (101) 
-Rl23 = —Rl23 — R213, Rl23 — R\23 — —R213, (102) 

M 1234 = M1234 = -M3413 = M4321. (103) 

We have chosen a rather general imaginary effective action at the required order which 
preserves gauge and chiral invariance, but we have included a larger number of terms than 
necessary to perform the matching process with the effective current. In fact, the matching 



22 



process could be done with solely the functions Qi 2 , P123, P123, and M 12 34. Instead, we have 
decided to include the additional term P 12 , in order to have the option of simplifying the 
action by a judicious choice of this extra function. For example, the extra function can be 
used to ensure that all functions remain finite at the coincidence limit, as will be explained 
later on. 

The calculation from the worldline formalism leads to the following contributions to the 
covariant current 

SW C = -i <T {ilPfPePJiSA, + i I 2 12 V a F^5A u + il^v a v a KD^mA v 
+ll 23 V^HV a V a H5A u + I^V^VaHVaHSAv + l\ 2 {D a UV ii V a mA v 

I 7 l2 ^^V a mA u + i I* 23 V a H^JA u + I^VpHVaHVaHSAv 
+I™ u V a HV f ,HV a mA u + I^DaHVaHV^HSAu) . (104) 

The coefficient functions are given in Appendix In order to express the current in this 
form, partial integration has been used to remove terms of the form DSA. In addition, 
indices that are contracted with the e tensor have been moved to the left, such that a term 
of the form VaV^ yields a sum of terms of the type VyV a and T^a. 

The contributions from the variation of Eq. (HOOD can be grouped in three levels, with 
the first level having only contributions from Q and P; the second level from the previous 
ones and R and R; the last level with all functions. Adding the contributions from the 
worldline method and the variation of Eq. fllOOp one obtains for the first level the following 
two equations 

P21 + {mi + m 2 )Q 12 = ll 2 , 
(mi + m 2 )P 12 - {ml - m 2 2 )Q 21 = if 2J (105) 



which have the solution 



^21 ; 2 



Qv* = 2 2 f^- ( 106 ) 

m\ — m 2 mi + m 2 

Besides, there arises the following restriction which is satisfied and can serve as a check for 
the corresponding terms in the effective current 

(mi + m 2 )ll x = -I 2 2 , Il 2 = -Il 2l l\ 2 = l\ 2 . (107) 
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The matching equations for the next level are 

- V 2 ((mi + m 2 )Q2i - Pgi) + Q12 + Qgi + (mi + m 3 )(-R 123 + R m ) = lf 23 , (108) 
V 2 ((mi + m 2 )(Q 12 + Q 2 i)) " (mi + m 3 )% 2 - 2Qu - 2Qgi + R312 = A 5 23 , (109) 
V 2 ((mi + m 2 )P 12 ) - (mi - m 2 )V 2 ((mi + m 2 )Q 21 ) 
-2Pi 2 + 2(mi - m 2 )g 2 i + (mi + m 3 )(Qi 3 + 2Q 3 i) 

-R 12 3(m 1 + m 3 ) + (m 1 -m 2 )(m 1 +m 3 )R3i2 = I 7 l23 , (110) 

and their complex conjugates. 

The first Eq. (11081) is of the form 

J3 23 

-R123 — -R312 = ; , (HI) 

mi + m 3 

and a set of solutions to Eqs. (11081) and (11091) is hence given by 

R 12 ^- l -{^-\ (112) 

2 V mi+m3 /l23 2 V mi+m V312 2 V mi+m3 /231 

P123 = %i - (mi - m 2 )i?i 23 . (113) 
The functions / and / are hereby defined as 

A23 = /i 3 23 + V 2 ((mi + m 2 )Q 2 i - P21) - Q12 - Q21, (H4) 

-fl23 = A 5 23 _ 

V 2 ((mi + m 2 )(Qi 2 + Q a i)) + 2Q i2 + 2Q 2i . (115) 

The last Eq. (jllOp leads to a constraint on the / functions that is given in Appendix O 

The function R possesses the required symmetries, and it reproduces the effective current 
correctly, but it is not necessarily finite in the coincidence limit, mi — > — m 3 . One way 
of solving this problem is to choose the function P appropriately which up to this point 
remained undetermined. Such a choice is e.g. given by 

P12 = — , Q12 = 0, (116) 
mi + m 2 

which leaves / as 

I 2 I 2 

J 123 = J 3 2 I =i (117) 

123 (mi - m 2 )(m 2 - m 3 ) (m 1 - m 3 ) (m 2 - m 3 ) 

With this choice, R is finite in the coincidence limit, as can be checked explicitly, and since 
/ is also finite, so is R. 
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For the last level, the following three equations hold 




■312 



2(mi + m 4 )Mi234 



I: 



■9 

12345 



(118) 



(V 3 - Vi) {Ranimi + m 3 )J + V 2 i?3i2 - 2i?3i 2 + 2^23 

-2(mi + m 4 )M4i23 
(Vi + V 2 ) f (mi + m 3 )^3i 2 ) + V 3 % 2 - 2^423 




(119) 



+2(m a + m 4 )Mi 2 3 4 ) 



-11 

1234- 



(120) 



One of these equations can be used to determine M, while the other two lead again to 
constraints on the I functions, the sum of the three equations has the especially simple form 



Since all previous functions in the effective action have been chosen finite in the coinci- 
dence limit, so is M1234. Eqs. (11 181) and (11201) show that M1234 is finite in the limit mi — * 1712, 
while Eq. (11191) shows that M1234 is finite in the limit ni\ — > — m 4 . This concludes the dis- 
cussion of the next to leading order contributions in two dimensions. 



We presented a worldline formalism to calculate the imaginary part of the covariant cur- 
rent of a general chiral model in the derivative expansion and our results are best summarized 



The resulting covariant current can be used to reproduce the imaginary part of the 
effective action by integration or matching. The advantage of this approach, compared to 
explicit formulas of the effective action, as given e.g. in ref. [191 ] . consists in the manifest chiral 
covariance. Chiral covariance reduces the number of possible contributions to the current 
enormously and makes even next to leading order calculations manageable as demonstrated 
in section IIVI in the case of two dimensions. Besides the chiral covariance, the use of the 
worldline formalism is essential in our approach. The evaluation of the worldline path 
integral requires neither performing Dirac algebra nor integrating over momentum space, in 
cont.asttothe^et— usedinref.fl. 



2(mi + m 4 )M 4 i 23 



(Vi + V 2 + V 3 )i?3l2 + A234 + A234 + I; 



11 

1234- 



(121) 



V. CONCLUSIONS 



by Eqs. ([62]) and ([63]) . 
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In principle, it is possible to use the presented formalism to determine the effective CP 
violation resulting from integrating out the fermions of the Standard model. For example, 
in next to leading order in four dimensions, an operator of the form VHVHTT could arise 
from the CP violation in the CKM matrix. Since the mass terms of the fermions are treated 
non-perturbatively in the derivative expansion, the resulting effective CP- violating operator 
is not necessarily proportional to the Jarlskog determinant. The discussion of this question 
is postponed to a forthcoming publication. 
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APPENDIX A: INTEGRALS USED IN THE CALCULATION 

In this section, the function g denotes the bosonic Green function 



Acknowledgments 



g(T,T!) = (y(T)y(n)), 



(Al) 



and 



g(T,T 1 ) = {y(T)y( n )) 



(A2) 



where the last expression does not contain a summation over the index A. 



1. Integrals in Two Dimensions 



In the calculation in two dimensions the following integrals have been used 




) 



(A3) 



The remaining occurring integrals can be expressed as 



1 1 1 

mf — m 2 



Jo T Jo Jo m 2 + m 3 



As = r% f dn r dT*e-^-^^^ g{T,r x ) V - J ' 7 
Jo 1 Jo Jo 

/CO jrp r~T fTl 
^ J dnj dr 2 e- T ^-^^^g(T,r 2 ) = -J. 



2 J 12 _ 

1 



m 2 + m 3 ' 



6 

321- 



2. Integrals in Four Dimensions 



In four dimensions the following integrals with three indices are used 



Jo ! Jo Jo 



J:. / f2 J o 1 J 2 ' 

m\ log(m 2 ) m\(m\ — mf) log(ra|) 



(m\ — m 2 )(ml — mf) (m\ — mf)(m\ — mf){m 2 — mf) 
ml log(m§) 



(m\ — mf)(m 2 — m\ 



2\ ' 



As = ^ [ dr 2 e-^-^^-^ g(T,n) 



m\ m\(m\ — m\mf) log(m 2 ) 



(mj — m 2 )(m 2 — mf) (m\ — m 2 ) 2 (m\ — mf) 
m\m\ log(m 2 ) m\m\ log(m|) 



[m\ — m 2 ) 2 (m 2 — mf) (mf — mf) 2 (m 2 — mf) ' 
As = [ % [ dn £ dr 2 e-^?"^- ■ (T> ^ = _ j; 
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The integrals with four indices can be expressed as 



Jo 1 Jo Jo Jo 

^3^123 



m 3 + m 4 ' 

"OO Jrp fT fT± PT2 



PCO Jrp pi pn PT2 

Jo 1 Jo Jo Jo 



V3J 9 



3^123 
7Tl 3 + m 4 ' 
-T 



J±L= f% \ dn r dr 2 r : dT3 e-^-^^-^ g(T,r 2 ) 
Jo 1 Jo Jo Jo 



^3^123 

m 3 + m 4 

fT 



pOO Jrp pi pri pT2 

Atu= / %~ 2 dn dr 2 dr 3 e- T ^-^^ + ^^(T,r 3 ) 
Jo 1 Jo Jo Jo 



- T 12 



Finally, the integrals with five indices read 

/■OO Jrp fT fT\ PT2 fTA 

4 5 345 = \ % dn dr 2 dr 3 / due -^-^^-^ 
Jo 1 Jo Jo Jo Jo 



V4J 11 



4^1234 

m 4 + m 5 ' 



4 6 345 = \ dn [ T1 dr 2 Pdr 3 rdne- T <^^-^g(T,n) 

Jo 1 Jo Jo Jo Jo 



V4J 12 



4 J 1234 



m 4 + m 5 

J\l^= / %- 2 \ dn dr 2 dr 3 dne~ T ^^+^g(T,T 2 ) 
Jo 1 Jo Jo Jo Jo 



V7 Tl3 

v 4^1234 

m 4 + m 5 

"OO jrp fT fT\ fT2 fT3 



/•CO jrp pi pTl pT2 pT3 

4 8 345 = / %\ dn dr 2 dr 3 dne- T ^^^-^ 9 (T,T 3 ) 
Jo 1 Jo Jo Jo Jo 



V7 714 

v 4^1234 

m 4 + M5 

"OO ^jrp r-T PTl PT2 PTi 



pCO jrp pi pTl pT2 PTZ 

4 9 345 = / % \ dn \ dr 2 / dr 3 / rfT4 e^-St^.« +1 -^ (r)T4) 
Jo 1 Jo Jo Jo Jo 



X7 I 12 

V 3^1234 

77l 3 + 7714 
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APPENDIX B: RESULTS IN FOUR DIMENSIONS 



In this section, we give the coefficient functions for the effective current and the effective 
density in four dimensions introduced in section IHI Dl 
The functions of the covariant current are given by 

2 morris — m\m 3 — m2,m 3 



1 123 



(mi + m 2 )(mi - m 3 )(m 2 - m 3 ) 

m\(mi(m2 — m 3 ) — 1m2m 3 ) \og[m\/m\^\ 
(mi + m2)(mi — m 3 ){m\ — m 2 v ){m\ — m 3 ) 
m2(m 2 (m 3 — mi) + 2mim 3 ) log [mj/mg] 



(Bl) 

(mi + 7YI2) {vrt2 — 7n 3 ) [m\ — m|) (m 2 , — m 3 ) ' 

^123 = — ^321) (B2) 



and 



^1234 = ^f 23 4 + ^234 log[m?] + A% U1 log^] + A^ 12 log[m 2 ] + A\ X23 log^] , (B3) 



with 



^1234 



mim27n 3 — m\m277i^ + m\m 3 m± — m 2 m 3 m4 
(mi — m 2 )(m 1 + 777.3) (mi — 777.4) (777.2 — ms)(m 2 + 7774) (m 3 — m 4 ) 



^ L mf(— — 777i?772 m3 + mim2m4 — 777,177737774 + 2777277737774) ^ 

1234 (mi — 777 2 ) (m! + 7773) (m! — 7774) (777^ — m|)(m 2 — m|) (777^ — m 2 ) 

The explicit functions occurring in the effective density are rather lengthy and hence we 
display them in terms of the integrals presented in the last section. 

5 i23 = J i23( m i + m s) - ^i 23 ( m i - m s) - Jllz{rn 2 - m 3 ), 
5 i234 = ^1234(^1 + m 4 ) - Jmt(mi - m 2 ) - Jll 3A (m 2 + m 3 ) + ^234(^4 + m 3 ), 
5 i234 = ^1234(^1 + m 4 ) - Jll^{mi + 777 2 ) - Ji2 3 34( m 2 - m 3 ) + J^ui^i + m 3 ), 

E>6 _ p>4 
-°1234 — ^4321 5 

5 1 7 2345 = ^12345 ( m l + m s) - ^345 ( m l + m l) + ^12345 ( m 2 + m 3 ) 

- J 12 8 345( m 3 + m 4 ) + ^I2 9 345( m 4 + (B5) 
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APPENDIX C: CO VARIANT CURRENT IN NEXT TO LEADING ORDER 



In this section we summarize the contributions to the covariant current in next to leading 
order. For the first two levels, they are given by 

I 2 

3 m 1 -m 2 + 2m 3 2 m 1 + m 3 4 

A23 = 2 2 V 2 (/ 2 i) + 8— 2 V 2 mim 2 J 12 

7?T.n — ?7iq - mi — mi 



+16 



2 2 
m 2 — m 3 



(mi -m 3 )(m 2 - m 3 ) ' 



r5 _ Q V 2( J 2l) o % .3 J 132 

123 — / \/ \ 123 ' 

mi — m 2 (mi — m 3 )(mi — m 2 J m\ — m 2 

,r . | / Vi/i 2 2 \ VJl 2 \ 3m 2 -m 3 2 

J = m 2 + m 3 ■ - + ■ V 2 I 12 

\mi+m 2 / 312 m 1 + m 2 l m 2 + m 3 



+4(m 2 + m 3 )(m 1 (m 1 —m 2 — m 3 ) + m 2 m 3 ) I I _Zi^2_ ) Vi J 12 

\\m 1 + m 2 J 3U m 1 + m 2 



-8(m 2 + m 3 ) 



'Vi(mim 2 J^ 2 )\ Vi(mim 2 J^ 2 ) 
m 1 + m 2 J 312 mi + m 2 



,/ . /V 2 ((m? - 4mim 2 + m 2 ) Jl) , \ T 4\\ /^,-,\ 

+4 m 2 - m 3 ^ ■ - V 2 mi - m 2 J x 4 2 ) , CI 

V m 2 + m 3 v 7 J 

and the relations 

r4 _ _t3 t6 _ _t5 t8 t7 fn^ 

J 123 ~~ J 321' J 123 ~~ J 321' J 123 ~~ J 321- V^^) 



30 



The contributions to the last level read 



7-9 _ _rll 

J 1234 ~~ J 432H 

rio . m i + m 4 .11 1 / m 4 + m 4 \ 5 

/l234 " _ ^^^ 123 + H V2 + (Vs) ^ 123 ) /l23 

+ (V 2 + V 3 + (V 3 ) 4123 ) ^^(4a - (V 2 4) 2 i 3 ) 

\ m 3 — m 4 m<2 — m?> — — 

+\ ( v 2 + V 3 + (-Vi + V 2 + V 3 ) 4123 ) (J 3 23 - V 2 I\ 2 ) 

Z \ 171% — m 4 - j 

+1 f-vi + v 2 + ^i±^i (v 3 ) 4123 ) (4i - (v 2 4) 321 ) 

2 \ m 3 — m 4 - / 

(-^213 ~~ (V 2 /i 1 2 ) 2 l 3 ) + — ——(-^432 _ (V 2 /i 2 ) 432 ) 



2(m 2 — m 3 )(m 3 — m 4 ) — — m 2 + m 4 

■(-^421 - (V 2 /i 1 2 )42l), 



2(m 3 — m 4 ) 

n 4 4(m 2 -m 2 )J 3 3 4(m 3 + m 4 )J 4 4 4J 4 5 24 

-'1234 - V 2-< - 



m 3 — m 4 (m 3 — m 4 )(mf — m 4 ) m\ — m\ m 3 — m 4 

mi + m 2 / m 2 + m 3 \ /Vi((m 2 - m^)(/ 1 2 2 - 8m 1 m 2 Jf 2 )) 

-v 3 



m 3 — m 4 \ m 3 + m 4 / \ (mf — m 2 )(mi + m 2 ) 
| ( m 3 + m 4 ) (Vi(7 2 2 - 8mim 2 J 4 2 )) 413 4(m 2 + m 3 )(m 1 -m A ) ^ jA 
(m 2 — ml)(mi — m 4 ) mi + m 2 )(m 3 — m 4 ) 1 

4(mi + m 2 ) / _ m 2 + m 3 \ (m 2 + m 2 ) J 4 3 - (m 2 - m 3 ) J, 



^23 



, -v 2 + — — — v, 

m 3 — m 4 \ m 3 + m 4 J mi — m 2 

4(m 2 + m 3 )(m 1 -m 4 )(m 3 + m 4 ) (ViJ 4 ) 312 4(mi + m 3 ) 5 

■ Vi ; V 3 J 



mi + m 2 m x — m 3 m 3 + m 4 

4(mi + m 2 )(m 2 + m 4 ) v ^ 6 ^ 



m 3 — m 4 



All functions / are finite in the coincidence limit and must fulfill the following constraints 
due to the behavior of the terms in the effective action under cyclic permutation and complex 
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conjugation: 

(m 1 +m 3 )/| 12 | ^ i| 

m 2 — m 3 321 (m 1 -m 2 )(m2-m 3 ) (mi — ■m 2 )(m 1 — m 3 ) 

(mi +m 3 )J 2 2 3 ijj = 

(mi — m 2 )(m 2 — m|) (m x — m 2 )(m 2 — m 3 ) 
3 .3 (mi + m 3 )/|i 

^123 + ^321 mi _ m2 + J 123 + ^321 

r2 r2 r2 

J 12 J 21 J 31 



J 12 _| _±\ ZJ?± 

■\ - m 2 (mi - m 2 )(m 2 - m 3 ) (mi - m 2 )(m 2 - m 3 ) 

(mi + m 3 )if 12 + (mi + m 3 )/| 12 - A 8 23 + 



mi — m 2 

+ ^1^±^4 = 0, (C4) 
(m x - m 2 )(m 2 + m 3 ) 



and 



A23 ~~ A23 ~~ 0' A23 ~~ A23 ~~ 0' A23 ~~ A23 ~~ 0- (C5) 
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